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Abstract 

The overall steady-state energy balance with two phases in a flow 
domain requires that the change in energy of the domain is equal 
to the difference between the total energy entering the domain and 
that leaving the domain. Prom the condition, the integral thermal 
flux across the surface is studied for a steady thermocapillary drop 
migration in a flow fleld with uniform temperature gradient at small 
and large Marangoni (Reynolds) numbers. The drop is assumed to 
have only a slight axisymmetric deformation from a sphere. It is 
identified that a conservative/nonconservative integral thermal flux 
across the surface in the steady thermocapillary drop migration at 
small/large Marangoni (Reynolds) numbers. The conservative flux 
confirms the assumption of quasi-steady state in the thermocapillary 
drop migration at small Marangoni (Reynolds) numbers. The noncon- 
servative flux may well result from the invalid assumption of quasi- 
steady state, which indicates that the thermocapillary drop migration 
at large Marangoni (Reynolds) numbers cannot reach steady state and 
is thus a unsteady process. 

Keywords Interfacial tension; Thermocapillary migration of drop; 
Quasi-steady state; Microgravity 
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1 Introduction 



The motion of a drop or bubble in the microgravity environment embedded in 
an immiscible mother liquid with a uniform temperature gradient is termed 
as thermocapillary migration of the drop or bubble, which is a very interest- 
ing topic for both fundamental theory and engineering application[l]. Young 
et al(1959) carried out an initial study in this area, called as YGB model[2], 
and gave an analytical prediction on its migration speed in the limit case of 
zero Reynolds(Re) and zero Marangoni(Ma) numbers, and a series of the- 
oretical analyses, numerical simulations and experimental investigations on 
this subject were carried out ever since. Subramanian(1981)[3] and Crespo 
et al(1996)[4] extended the YGB results to small Ma numbers and obtained 
analytical results in series expansion of Ma numbers. With consideration of 
thermal boundary layer, the analytical results for migration speed of a bub- 
ble at large Ma (Re) numbers[5,6] agree well with the corresponding results 
of steady state numerical simulations[7,8,4] and experimental studies[9]. 

Although the thermocapillary bubble migration processes are understood 
very well, the behavior of thermocapillary drop migration appears rather 
complicated due to the transfer of momentum and energy though the in- 
terface of two-phase fluids. For the migration of a drop, on the one hand, 
the experimental result of the migration speed at small drops 11 ± 1.5//m in 
diameter obtained by Braun et al(1993)[10] agrees with the YGB model. On 
the other hand, another experiment for larger drops with diameters rang- 
ing from 0.69 to 2.38 mm performed by Wozniak(1991)[ll] shows that the 
migration velocities arc smaller than those given by the YGB linear pre- 
diction. Afterward more attention has been paid to thermocapillary drop 
migration for large Ma (Re) numbers. Hadland et al(1999)[9] carried out 
experiments based on Fluorinert liquid FC-75 and lOcst sihcone oil as the 
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drop phase and continuous phase during the NASA Space Shuttle mission 
with the maximal Ma(Re) number up to 3300(49.1). It was shown that the 
drop migration speed nondimensionalized by the YGB velocity decreased 
as the Ma number increased and the global migration process exhibited an 
unsteady nature. To further observe the variation trend of drop migration 
velocity with increasing Ma number, Xie et al(2005)[12] adopted Fluorinert 
liquid FC-75 and 5cst silicone oil as the drop and the continuous phases, re- 
spectively, and performed experiments in the Chinese spacecraft ShenZhou-4. 
The experimental investigation was completed for several ranges of large Ma 
(Re) numbers, where the drop didn't reach the steady state in the migration 
process, and the maximal Ma(Re) number reached 5525(302.6). It was also 
observed that the non-dimensional drop migration velocity decreased as the 
Ma numbers increased. However, from the theoretical analysis for the large 
Ma (Re) numbers, it was reported [13] that the migration speed of a drop 
increased with increasing Ma number, which is in qualitative agreement with 
the corresponding numerical simulation[14]. Both the theoretical analysis 
and numerical simulation are based on the assumptions of quasi-steady state 
and non-deformation of the drop. The above qualitative difference between 
experimental observations and theoretical/numerical results may result from 
the quasi-steady state or non-deformation assumptions of the drop in the 
model. Moreover, Herrmann et alfl5j and Wu et al[16j adopted respec- 
tively the numerical methods to investigate the thermocapillary motion of 
deformable and non-deformable drops and indicated that the assumption of 
quasi-steady state was not valid for large Ma numbers. Therefore, the ther- 
mocapillary drop migration at large Ma (Re) numbers is still a topic to be 
further studied with emphasis laid on its physical mechanism. 

To address the discrepancies between experimental and theroretical/numerical 
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results, in this paper, our effects are focused on the assumption of quasi- 
steady state in the process of thermocapillary drop migration. The drop 
may have only a slight axisymmetric deformation from a sphere. By using 
the asymptotic expansion method, we investigate the continuity of integral 
thermal flux across the surface based on the overall steady-state energy bal- 
ance in the flow domain, and analyze the existence of quasi-steady migration 
of the drop at zero, small and large Ma (Re) numbers. 



2 Models and quasi-steady state assumption 



Consider the thermocapillary migration of a spherical drop of radius Ro, 
density jp, dynamic viscosity afi, thermal conductivity /3k, and thermal 
diffusivity Xk in a continuous phase fluid of infinite extent with density p, 
dynamic viscosity /i, thermal conductivity k, and thermal diffusivity k under 
a uniform temperature gradient G. The change rate of the interfacial tension 
between the drop and the continuous phase fluid with temperature is denoted 
by (Tt- Axisymmetric energy equations for the continuous phase and for the 
fluid in the drop in a laboratory coordinate system denoted by a bar are 
written as follows 



where v and T are velocity and temperature, and a prime denotes quantities 
inside the drop. The solutions of Eq. (1) have to satisfy the boundary 
conditions at infinity 



§ + vVT = kAT, 
^ + v'Vf ' = XkAT', 



(1) 



^ To + Gz, 



(2) 
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where Tq is the undisturbed temperature of the continuous phase and the 

boundary conditions at the interface {fb, Zb) of the two fluids 

f{n,zb,t)^f{n,zb,t), 

^{fb,zb,t) = /3^{n,zb,t), 

where n is a unit vector normal to the interface. In what follows, the undis- 
turbed temperature Tq is reduced for simplicity. 

In general, the surface tension decreases with the increasing of the local 
temperature. For a temperature field with its gradient in the z direction, the 
generated surface tension force is a net force along the surface and the droplet 
starts to move towards the warm side under the action of net force. When the 
net force acting on the drop at the flow direction is zero, the thermocapillary 
drop migration reaches a stable process. However, due to the variation of 
physical parameters with the ambient temperature, the migration process 
may not reach stable state. Only when the migration is sufficiently slow 
that the order of relevant time scale for the transport process to generate 
stable velocity and temperature fields is smaller than that for the drop to 
move an appreciable distance, the assumption of the quasi-steady state is 
valid. It means that after experiencing an initial unstable migration process, 
the drop migration may reach a steady state at the time to and the position 
Tq — Zok, i.e., migrating with a constant drop migration speed V^. Using 
the coordinate transformation from the laboratory coordinate system to a 
coordinate system moving with the drop velocity V^x, 

f^r + ro + V^{t-to)k, v(f , = v(r) + Kok, f{f,t)^T{r) + G[zo + Vo^{t-to)], 

v'(r, t) = v'(r) + Kok, f'{r, t) = T'(r) + G[zo + V^it - to)], 

(4) 
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the problem (1) can be formulated as 

GKo + vVT^/^AT, 

(5) 

GV^ + v'VT' = AkAT'. 
The details of the transformation are given in Appendix 1. By taking the 
radius of the drop i?oi the velocity Vq = —aTGRo/ fi and GRq as the reference 
quantities to make the coordinates, velocity and temperature dimensionless, 
energy equations (5) combined with the continuous equations can be written 
in the following dimensionless form in a spherical coordinate system (r, 9) 

V^ + V- (vT) = ^ AT, (6) 

+ V ■ {VT') = ^AT', (7) 

where v = {u, v) and Marangoni number is defined as 

Ma = (8) 

K 

By using the transformation (4), the boundary conditions (2) and (3) can be 
respectively written in the form of dimensionaless as follows 

T — >■ r cos^, as r — >■ oo (9) 

at places far away from the drop and 

T{ro,e) = T'{ro,9), (10) 
dT dT' 

= (11) 

at the interface of the two fluids. Thus, once the drop migration reaches a 
steady state, the above problem in the laboratory coordinate system can be 
described by the steady energy equations (6) (7) with the boundary conditions 
(9) (10) (11) in the coordinate system moving with the drop velocity. This 
implies the overall steady-state energy balance with two phases in the flow 
domain in the co-moving frame of reference. 



3 Integral thermal flux across the drop sur- 
face under the quasi-steady state assump- 
tion 

In general, for a two-phase flow, as is the case in the present problem, if the 
quasi-steady state assumption is valid, the solutions of the problem not only 
satisfy the differential energy equations with boundary conditions, but also 
the overall steady-state energy with two phases in the flow domain under 
integral boundary conditions keeps balance. However, if the quasi-steady 
state assumption is invalid, the overall steady-state energy with two phases 
in the flow domain under integral boundary conditions is not balanced. This 
means the solutions of the problem also cannot satisfy the differential energy 
equations with boundary conditions. Thus, to confirm whether the thermo- 
capillary drop migrations at different Ma (Re) numbers are always in the 
quasi-steady state processes, we may analyze the overall steady-state energy 
of two phases in the flow domain in the co-moving frame of reference under 
integral boundary conditions. 

For the thermocapillary drop migration, to get the overall steady-state 
energy transport of two phases in the flow domain in the co-moving frame of 
reference, we have to integrate Eq. (6) and Eq. (7) in the continuous phase 
domain (r e [ro,roo],^ £ [O?'^]) within the drop region (r e [0,ro],^ G 
[0,7r]) as 

poo PTT poo p7V 1 poo p7V 

/ / V^dV+ / V-(vT)dV=— / ATdV, (12) 

Jro Jo Jro Jo Ma Jrg Jo 

rro r-TT rro rir X rro rir 

/ / V^dV+ / V-{v'T')dV^— / AT'dV, (13) 
Jo Jo Jo Jo Ma Jo Jo 
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and Eq. (10) and Eq. (11) at the drop surface (r — ro,0 E [0, tt]) as 

rT{ro,e)dS^ rT'{ro,e)dS, (14) 
Jo Jo 

where dV — r'^sm.OdrdO and ^5" = r^sinOdO. And then transforming the 
volume integration of Eq. (12) and Eq. (13) in the flow domains to the 

surface integration over the droplet surface and the surface at infinity in 
terms of the Gaussian formula, we have 



3 

and 



{rl-^ I rlsm0d0)+f{uT)\r^dS-f{uT)US = —{f ^V^dS-j ^\r,dS) 

(16) 



-fl .„3,in« + /(«'r')|..dS = — /-|.,dS. (17) 

Using the zero normal velocity boundary condition at the interface, we can 
derive 

^U^sinede ^ rlJ,-§\r^8mede-MarlJ^{uT)\r^smede ^^g^ 
-^(r^-l/oVgsin^de) 

and 



/ 

^0 



T 8T' V Mn f-^ 

l^rlsmede ^ r^smede, (19) 

3A Jo 

where the outer normal vector at infinity is the radial coordinate axis. Thus, 
Eq. (18) and Eq. (19) display integral thermal fluxes across the drop surface 
obtained from the overall energy transport. We assume that the drop has 
only a slight axisymmetric deformation from a sphere 

ro = l + /(e), /«1. (20) 
For this case, Eq. (18) and Eq. (19) may be written as 

J^^Ul^mede^ rl!,^ §\r^smede - MarlJ^iuT)\r^smede ^^^^ 

-'-^[rl - 1 - i /oV sin 0d0 + Oif)] 



dn 



rlsmede^-^^[2 + 3 f smOde + 0{f)]. (22) 
3A Jo 



and 

rdr 2 . V^Ma 
Jo 

To next-to-leading (first) order in /, Eq. (21) and Eq. (22) may finally be 
written as 

_2V^(^3^_l_|j;/si„^d^) 

and 

/ ^|,„r2sin^de = ^^(2 + 3/ fsmOdO). (24) 
Jo on iX Jo 

Prom Eq. (14) and Eq. (15), we have integral boundary conditions across 
the drop surface 

r T\ryoSm9d9 ^ ^ T'\r,r'^ sin 9d9, (25) 
Jo Jo 

—\r,Tl Sin 9d9 -f^l krl sin 9d9. (26) 

Thus, for a quasi-steady state thermocapillary migration of the drop with the 
slight deformation from the sphere, the overall steady-state energy balance 
of two phases in the flow domain in the co-moving frame of reference requires 
that the integral thermal fluxes (23) and (24) at the drop surface obtained 
from the overall energy transport are self-consistent with the integral bound- 
ary condition (26). In the following, we will investigate the self-consistency 
for the different Ma(Re) numbers. 

3.1 Conservative integral thermal flux across the drop 
surface at zero Ma (Re) numbers 

In the case of zero Re {Re = ^) and zero Ma numbers, i.e. the YGB 
model, scaled velocity and temperature fields of the continuous phase and 
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within the drop in Eq. (6) and Eq. (7) may be described[21 [IZ] as 



u 



V 



V^sme{l + ^), (27) 



and 



T =(r + i4^)cos^, 



u' = |KoCose(l -r2) 



v' = 

^ 2 



V;oSine(l-2r2), (28) 



T' = 2+^^ COS 9. 

Since Ma = 0, using the temperature field in (27), we can derive the following 
equality from Eq. (23) and Eq. (24) 

/ T^lro^nSin^rfe = /3 / —L.rlsindde = (29) 
Jo on Jo on 

Thus, under the quasi-steady state assumption, the integral thermal fiux 
across the drop surface at zero Re and zero Ma numbers is conservative, 
which corresponds to the integral thermal fiux boundary condition (26). This 
implies the overall steady-state energy balance of two phases in the fiow do- 
main in the co-moving frame of reference. The quasi-steady state assumption 
is valid. 



3.2 Conservative integral thermal flux across the drop 
surface at small Ma (Re) numbers 

For small Re number, the velocity fields in Eq. (6) and Eq. (7) may be 
described by the creeping fiow. The general solutions of the scaled fiow field 
in the continuous phase and within the drop are given by [ISl [H] as 

U =-V;oCOS^(l-^)-(l-^)E~3^nr-"+lPn-l(cOS^), ^^^^ 

V = KoSin^(l+ 2^3) + E^=3^n(7^ - 7^)C„-^/'(cos^^)/sin^ 

11 



and 



u 



I _ 3 



fV;, cos 0(1 - r2) + - 1) E:L3^nr"P„_i(cos( 



v' 



-|\4oSin0(l - 2r2) - ^^^3 D„[nr"-2 - + 2)r"]C-^/2(cos 0)/ sin 0, 

(31) 

where C^^/^(cos6') is the Gegenbauer polynomial of order n and degree — 
Pn_i(cos6') is the Legendre polynomial of order n. Dn is given as 

And scaled temperature fields in the continuous phase and within the drop 
at the small Ma numbers are given in [19] as 

+ 3A(2W^(2+3a) [^ + ^ + ^2(cOS0)(f + f| + ^)]e + 0(6^), 

+ A(2W(2T3^['^i + '^2^' - t^' + ^2(cos0)(5^r2 + fr^)]e + ©(e^) 

(33) 

where 5i = 2[A(1 - /3) - /3(2 + /?)], ^2 = -|(1 - /?), ^3 = -A(4 - /3), 
^4 = 7(3T2^[7A(8 + 5/3 - 4/32) _ ^g^j^ ^/ ^ X[6A(1 - /3) - (8/3^ + 20/3 + 17)], 
5^ = 1(2/3 + 13), 5^ = - 21(3+2/3) [7A(7-/3) + 9(3+4/3)] and the small parameter 
is e = Ma. Then, using the temperature field in (33), we simplify Eq. (23) 
and Eq. (24) to 

JO dn\ro' 0^^^^^^^ ~ 3A(2+/3)2(2+3a)'^"l ^ rg^'"' ^ 3{2+/3)'^-^ rl^> 

+V^Jofsmedee + 0{e^) 
^ - 3A(2+;A2+3.) e + l0jV^e + V^ J- f sin OdOe + 0(6^ 

(34) 

/ -T^Lrl sin 9 d9 = ^{2 + 3 fsmede)e, (35) 
JO on OA Jo 
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and 



where = (2+/3)^2+3a) migration velocity of the droplet at zero Re 

and zero Ma numbers. From Eq. (34) and Eq. (35), we have 

Iomro-mroyosin9de= Uf-'^)iV^-V^)e ^^^^ 

+ (l-f)\4o/o"/sin» + 0(e2). 

Since e 0(\4o V^) and / < 1, we have 

/ T^LrUinede ^ (3 —lynSmOde. (37) 
Jo on Jo on 

Thus, under the quasi-steady state assumption, the integral thermal flux 
across the drop surface at small Re and small Ma numbers is conservative, 
which corresponds to the thermal flux boundary condition (26). This implies 
the overall steady-state energy balance of two phases in the flow domain in 
the co-moving frame of reference. The quasi-steady state assumption is valid. 

3.3 Nonconservative integral thermal flux across the 
drop surface at large Ma (Re) numbers 

For large Re number, the velocity fields in Eq. (6) and Eq. (7) can be 
described by potential flows and boundary layer flows [211]. The scaled inviscid 
velocity field in the continuous phase and Hill's spherical vortex within the 
drop can be respectively written as 

u = — yooCos6'(l — \), 

V ,3;, ^gg^ 

V = VooSm9il + ^) 

and 

u' = ^ cos 6^(1 - r^), 

^ ^ ' (39) 
v' = -^sin^(l - 2r2). 

Since Eq. (9) only gives the primary approximation of the temperature field 
at infinity, we have to obtain an asymptotic expansion of T for the integration 
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of Eq. (23). To determine the asymptotic behavior of T at r ^ 1, the 
analytical result of outer temperature field in the continuous phase at the 



small parameter e = 1/ y/Vo^Ma is given in [6] as 

T = rcos6'+ / (f sin6' — M cos6' — l)/'u|^c/f + o(l), (40) 

J oo 

where \E'[= |sin^^(r^ — 1/r)] is the streamfunction of the continuous phase 
and the symbol " +" before the integral is determined to preserve the monotonously 
increasing trend of T(r, 0) with r(> 1) in the continuous phase. Using Eq. 
(38), it can be derived as: 

r 1 1-3 sin^ ^ , , 
T = r cos ^ + ^^^kc^r + o(l). (41) 



T = rcose+l \ :1 . .. Wdf + ojl), (42) 



Replacing 6' by \E' in Eq. (41), we have 

1 3^/(f2 - l/f) - 1 
r-' - 1 ±^l-2^/(f2_l/f) 

where the symbol " ±" in the integral depends on the value of 6 (the symbol 
" + "/" - " corresponds to G [0, 7r/2)/[7r/2, vr)). At r > 1, the result (42) 
can be expressed as 



T ^rcos^ + /,°°^ \^df + o{l) 



r cos 6 — 1^ cos 6 + o{l) 



(43) 



where ~ | sin^ Or"^. 



Using the temperature field at the infinity in (43), we can simplify Eq. 
(23) and Eq. (24) and derive 



/j^ f \ry, sinede = -^(1 - -^) + ^ £ f sin ede + o(^ 



and 

f^dT'. o „.„ 1 



(44) 



/ —\r,rlsmede = -—{2 + 3 fsinede). (45) 
Jo on 3Ae^ Jo 
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f7\ 

Jo 



L„rn sin 6'd6' > / — L„ sin 6'(i6'. (47) 
Jo on 



Prom Eq. (44) and Eq. (45), we have 

/o^(/5£l.o - f IJro^sin^d^ = 3^,(1 + f ) + (f - l)i/o'^/sin^d^ 

= 1(1 + f )KoMa + (f - l)V^Ma f sin 9d9. 

(46) 

Since both P and A are positive and / <^ 1, we have 

dT'. , . „ .„ dT. 
dn 

So, if the overall steady-state energy with two phases in the flow domain 
under integral boundary conditions is balanced, Eq. (47) should be reduced 
to Eq. (26), which seems impossible. It is termed as a nonconservative in- 
tegral thermal flux across the surface for the steady thermocapillary drop 
migration at large Ma (Re) numbers. This imphes the overall steady-state 
energy unbalance of two phases in the flow domain in the co-moving frame of 
reference and indicates that the thermocapillary drop migration at large Ma 
(Re) numbers cannot reach steady state. Thus, it is clear that the invalid 
assumption of quasi-steady state for the thermocapillary drop migration pro- 
cess is a reasonable explanation for the nonconservative integral thermal flux 
across the drop surface. 

To analyze the thermal flux near the boundary, we write the integrals of 
Eq. (44) and Eq. (45) in the discretization scheme as follows 

and 

Jo 9^ ^"^^^^^ ^'^^ ^ °' ^^^^ 

where 6i G [0,7r] and A^^ = n/N. Since rpsin^^j > 0, we reach a conclusion 

that there must be some interface points 9i e [0, tt] where the following 
equation holds 

dT dT' 

-d^M)<^<^M). (50) 
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or some interface points 6i and 6j G [0, tt] where the following equations hold 

Physically, this means that near these points 9i the thermal energy is trans- 
ferred from the interface to outside (the surrounding fluid) as well as from 
the interface to inside (the droplet) or near these points 0i/9j the transfer- 
ence of thermal energy from outside/the interface to the interface/outside is 
weaker /stronger than that from the interface/inside to inside/the interface. 
On the one hand, if Eq. (50) can satisfy the integral thermal flux boundary 
condition in Eq. (26), thermal sources inside the interface will be introduced 
to balance the transference of thermal energy. On the other hand, if Eq. 
(51) can satisfy the integral thermal flux boundary condition in Eq. (26), 
thermal sinks inside the interface or thermal sources in the droplet will be 
introduced to decrease the transference of thermal energy from the interface 
to outside or increase the transference of thermal energy from inside to the 
interface. Since there is absolutely no thermal sources or sinks inside the 
interface or thermal sources in the droplet, the above transport processes of 
thermal energy near the interface seem impossible. It means that the thermal 
flux across the drop surface is nonconservative. 

4 Conclusions 

In summary, from the condition of overall steady-state energy balance with 
two phases in a flow domain, we have identified a conservative/nonconservative 
integral thermal flux across the surface for a steady thermocapillary migra- 
tion of a drop with a slight axisymmetric deformation from a sphere in a 
uniform temperature gradient at small/large Ma (Re) numbers. The conser- 
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vative integral thermal flux conflrms the assumption of quasi-steady state in 
thcrmocapillary drop migration at small Ma (Rc) numbers. The nonconser- 
vative integral thermal flux may well result from the invalid assumption of 
quasi-steady state, which indicates that the thcrmocapillary drop migration 
at large Ma (Re) numbers cannot reach steady state and is thus a unsteady 
process. 
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and IMECH/SCCAS SHENTENG 1800/7000 research computing facihties 
for assisting in the computation. 



17 



Appendix 

Based on the transformation between two cylindrical coordinate systems 
(f , z) and (r, z) of Eq. (4) , we have 



^ It 9r 1*^ ~^ 9z I 9r 1*^ ~^ 92 ^ 1 1 ' (52) 

Alt = mri) + iirmt = + iirB]\t = Al*. 

And for energy equation Eq. (1) of the continuous phase fluid, we can also 
write its unsteady, convection and conductivity terms as follows 

df I — dT\_ 4_ r<\/ — dT\ dr\_ , dTl a2|__|_9T| dt\_ , /-«t/ 

Qtlr — dt \^ ' ^^'^ ~ dr dt\^ ' dz I* Stir dt li"Stlr "'"^''00 

vVf|i ={v + V^k)V{T + GV,,t)\t = {v + V^k)VT\t (53) 

= vVT + Kof, 
AT I, = A(r + GKoi) It = ATI* = AT, 

where |^|r = f^lf = 0, |f If = |f |j = -Kx> and ^|r = 0. Then, substituting 
Eq. (53) into Eq. (1) in the laboratory coordinate system, we obtain Eq. (5) 
in a coordinate system with its origin fixed at the drop center 

GV^ + vVT = kAT. (54) 

Similarly, we can also transform the energy equation within the drop as 
above. 
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